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ABSTRACT 

Very early optical and near infrared (IR) emission was discovered accompanying the long gamma-ray 
burst (GRB) 041219a. We show that the optical/IR flash tracking the gamma-ray lightcurve during 
the prompt emission could be understood as emission from neutron-rich internal shocks, as has been 
suggested by Fan & Wei. The early if s -band afterglow lightcurve after the prompt phase could be 
well- modeled as the superposition of a reverse shock and a forward shock component. The radio data 
also support the reverse shock interpretation. 

Subject headings: Gamma Rays: bursts— ISM: jets and outflows-radiation mechanisms: nonthermal 



1. INTRODUCTION 

GRB 041219a was detected both by the IBIS (Imager 
on Board the INTEGRAL Satellite) detector of the IN- 
TEGRAL satellite (Gotz et al. 2004) and by the Swift 
Burst Alert Telescope (BAT) (Barthelmy et al. 2004). 
This burst distinguishes itself from other bursts in sev- 
eral aspects. (1) It is very bright. The 15-350 keV 
fluence measured by the Swift BAT was approximately 
1.55 x 10~ 4 ergs/cm 2 (Barthelmy et al. 2004), placing it 
in the top few percent among the 1637 GRB events listed 
in the comprehensive fourth Burst and Transient Source 
Experiment (BATSE) catalogue (Paciesas et al. 1999). 
(2) The duration of prompt 7— ray emission (T 90 ) is ap- 
proximately 520 seconds, making it one of the longest 
bursts ever detected. (3) The prompt optical and in- 
frared (IR) emission was detected to accompany the 
prompt 7— ray emission (Vestrand et al. 2005; Blake et 
al. 2005). This is the second case since GRB 990123 (e.g., 
Akerlof et al. 1999). (4) The location of this burst is near 
the Galactic plane and in a direction with high optical 
extinction (Galactic coordinates 1 = 120°, b — +0.1°), 
so that the R-band extinction is very large (~ 4.9 mag- 
nitudes or even larger). For this reason, no late optical 
afterglow has been detected. Fortunately, in the IR band, 
the afterglow (including that in the very early phase) has 
been well detected (Blake et al. 2005). The redshift is 
unknown. In this Letter, we assume z = 1. 

The very early long- wavelength observation is very im- 
portant to diagnose the outflow composition (Fan et al. 
2004; Zhang & Kobayashi 2005; Fan et al. 2005) \ since 
the late afterglow taking place hours after the burst trig- 
ger is powered by the forward shock (FS), so that essen- 
tially all the initial information of the ejecta is lost. In 
this Letter, we apply our previous analyses to discuss the 
very early long-wavelength observation of GRB 041219a. 

2. THE PROMPT OPTICAL AND NEAR-IR FLASH 

1 In principle, the ejecta of GRBs may be Poynting flux domi- 
nated (see Lyutikov & Blandford 2003 and the references therein) 
or neutron-rich (e.g., Derishev et al. 1999; Beloborodov 2003; 
Pruet et al. 2003) or both (e.g., Vlahakis et al. 2003). 



There has been some interest in discussing/searching 
for prompt long-wavelength radiation accompanying 
prompt 7— rays emission even in the pre-afterglow era 
(e.g., Katz 1994; Schaefer et al. 1994; Wei & Cheng 
1997). In the afterglow era, more theoretical attention 
was paid on the topic. In the standard internal shocks 
model, accompanying the prompt 7— ray emission, long 
wavelength flashes are expected (Meszaros & Ress 1997, 
1999; Fan & Wei 2004a). If there are a large amount 
of neutrons contained in the GRB outflow, the decayed 
neutron shells would provide more collisions at a larger 
distance from the central engine. If a burst is long 
enough, the neutron-rich internal shocks would give rise 
to detectable long wave-length flashes during the prompt 
7— ray emission phase (Fan & Wei 2004b). 

In the standard internal shock model, the 
synchrotron-self-absorption frequency could be 
estimated as (e.g., Li & Song 2004) z/ a ~ 

7 x 10 16 Hz Ll(l Ko rTo, RryXB'V 7 ^), where 



L syn is the synchrotron radiation luminosity, B' is the 
comoving-frame magnetic field strength in the internal 
shock phase, i?i n t ~ 2r i 2 n c<5i/(l + z) is the typical internal 
shock radius, r m is the bulk LF of the merged two 
shells, St is the observed variable timescale of the prompt 
7— ray lightcurve, and the convention Q x = Q/10 x has 
been adopted in cgs units here and throughout the 
text. For typical GRB parameters, one can see that 
v & is usually well above the optical band. This tends 
to suppress the optical flux. Also the F v spectrum is 
expected to have a power law index 5/2, inconsistent 
with the prompt IR data of GRB 041219a (see Figure. 2 
of Blake et al. 2005). Although with proper adjustment 
of parameters, the proton-dominated internal shocks 
models may be able to match the observation, in this 
Letter we focus an alternative interpretation, i.e. the 
neutron-rich internal shock model (Fan & Wei 2004b). 

There are good reasons to assume that a large amount 
of neutrons (comparable to the amount of protons) ex- 
ist in the GRB ejecta (e.g., Derishev et al. 1999; Be- 
loborodov et al. 2003; Pruet et al. 2003). In the neutron- 
rich internal shock model, the LFs of the proton shells 



2 



are variable, so are the LFs of accompanying neutron 
shells. For the slow neutrons (with r niS = 50) coupled 
with the slow proton shells, they do not interact with 
other materials before decaying (the typical 0— decay ra- 
dius is ~ 900 s oT n . s = 1.3 x 10 15 cm, where c is the 
speed of light). However, if the GRB lasts long enough, 
the slow neutron shells ejected at earlier times would 
be swept successively by the faster proton shells ejected 
at later times. This happens in the distance range of 
~ 10 13 — several x 10 15 cm, in which the slow neutron 
shells decay continuously. The proton shells interact with 
the (3— decay products of the slow neutron shells and 
power detectable long wavelength prompt emission, as 
shown in Fan & Wei (2004b). 

A detailed treatment of the process has been presented 
in §2 of Fan & Wei (2004b). A novel effect taken into ac- 
count here is the inverse Compton (IC) cooling of the 
electrons because of the space-time overlapping between 
the proton shell— neutron shell interaction region and the 
prompt MeV 7— ray photon flow (e.g. Beloborodov 2005; 
Fan et al. 2005). The inverse Compton parameter is cal- 
culated as Y = P IC /P syn , where P IC is the energy loss 
rate of one electron via IC scattering with the prompt 
7— ray emission, P syn = (4/3)ctt7c/3c ^B c is the energy 
loss rate of one electron via synchrotron radiation, 7 C 
is the random Lorentz factor of the emitting electron, 
and Ub is the magnetic energy density generated dur- 
ing the interaction. Since the Klein-Nishina correction is 
important for the IC process discussed here, we strictly 
use Eqs. (2.47-2.51) and Eq. (2.56) in Blumenthal & 
Gould (1970) to calculate P IC , The prompt 7— rays pho- 
ton number distribution is taken as n ey = n e b (e 7 /e 7 )~ 3 / 2 
for e^/50 < e 7 < and n ty = n e b(e 7 /e^y 2 - 25 for 
e 7 < e 7 < 50e 7 , where e 7 ~ 250 keV is the peak en- 
ergy of the observed spectrum 2 , n e h is constrained by the 

condition f n e e 1 de 1 = U 1 , where \J 1 w L 1 /AnR 2 ^ 2 c is 
the initial 7— ray photon energy density, 7 is the evolving 
bulk LF of the proton shell during the sweeping process 3 , 
and L-y is the prompt 7— ray luminosity. 

We calculated the synchrotron lightcurves for a proton 
shell interacting with the decay products of a series of 
slow neutron shells. The result is shown in Fig. ^ Fol- 
lowing parameters have been taken: M p = 3.7 x 10 26 g 
is the rest mass of one proton shell; M® — 10 26 g is 
the initial rest mass of one neutron shell; r m = 600; 
r n .s = 50 (With these parameters, we get 7 ~ 300 at 
the end of the interaction. This will be regarded as 
the intial Lorentz factor in the afterglow phase, and we 
will rewrite it as 77 in ^31 This matches the one found 
in our modeling the early K s band afterglow. See SJ2J); 
St = 10~ 2 (1 + z)s; D L = 2.2 x 10 28 cm is the luminosity 
distance; L 7 = 10 52 ergs s _1 ; and e c = 0.3 and eb = 0.1 
are the fraction of shock energy given to the shocked elec- 
trons and magnetic field, respectively. The lightcurves 
reach the peak as the observer frequency v b s crosses the 
synchrotron-self- absorption frequency 24,, and another 
break during the decaying phase marks the epoch as v h s 
crosses the typical synchrotron radiation frequency. The 

2 Given a same e 7 , the calculated results are not very sensitive 
to the photon spectral indices assumed. 

3 Notice that we have adopted the same symbols as in Fan & 
Wei (2004b) to keep consistency between the two papers. 
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Fig. 1. — Long wavelength (the observer frequency has been 
marked in the figure) emission powered by a proton shell inter- 
acting with the decay products of a series of slow neutron shells 
as a function of time (Similar to Fig. 2 of Fan & Wei (2004b)). 
Following parameters are taken: M p = 3.7 X 10 26 g, Mj = 10 26 g, 
T m = 600, r n , s = 50, St = 10^ 2 (1 + z)s, D L = 2.2 x 10 28 cm, 
£ C = 0.3, e B = 0.1, L-y = 10 52 ergs s" 1 and = 250keV. The 
real lightcurve should include the contributions from many proton 
shells. 



v a is much smaller than that in the standard internal 
shock model because of the following reasons. Firstly, 
we have R ^> i?i n t- Secondly, the forward shock up- 
stream proton number density is much smaller than that 
in the standard internal shocks (see eq.[l] of Fan & Wei 
(2004b) for detail), which results in a smaller B' and 
weaker synchrotron emission. Finally, with the IC cool- 
ing effect taken into account, and since synchrotron emis- 
sion is weaker, the Y parameter is large, i.e. about tens to 
hundreds with the typical parameters. The synchrotron 
radiation luminosity is further lowered by a factor l/Y. 

The detected flux should be the integrated emission 
powered by a series of proton shells interacting with 
the decay trail of slow neutron shells. For j/ b s — 
(1.4, 1.8, 2.2, 4.3, 10.0) x 10 14 Hz, the predicted flux are 
(18.8 ,21.2 , 23.3, 26.3, 27.2) mJy, respectively. Such 
strong emission is detectable with the current telescopes, 
and is consistent with the optical observation of GRB 
041219a when the extinction correction is taken into ac- 
count (Vestrand et al. 2005). The averaged spectrum 
is flat, which is roughly consistent with the earliest IR 
band observation of GRB 041219a (Blake et al. 2005). 
As mentioned in Fan & Wei (2004b), the predicted long- 
wavelength emission is expected to be correlated with the 
prompt 7— ray emission but has a ~ 10(1 + z)T~l 17 s 
lag. This is also consistent with the tracking behavior of 
the prompt optical flashes (Vestrand et al. 2005). Notice 
that it is not our intention to fit the prompt optical and 
near IR emission lightcurves closely (which requires addi- 
tional complicated assumptions about the central engine 
behavior). The main purpose of the current discussion 
is to indicate that the neutron-rich internal shock model 
suffers less constraints than the proton-dominated inter- 
nal shock model, and can better account for the observed 
prompt optical/IR spectrum. 

3. THE VERY EARLY NEAR-IR AFTERGLOW 

After the internal shock phase, as the fireball is deceler- 
ated by the circumburst medium, usually a pair of shocks 
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develop (Meszaros & Rees 1997; Sari & Piran 1999). The 
early optical afterglow lightcurve is usually composed of 
the contributions from both the forward and the reverse 
shocks. Zhang et al. (2003) pointed out that depending 
on parameters, there are two types of early optical/IR 
lightcurves for a fireball interacting with a constant den- 
sity medium (ISM), i.e. Type I (rebrightening) in which 
distinct reverse shock and the forward shock peaks are 
detectable, and Type II (flattening) in which the forward 
shock peak is buried beneath the reverse shock peak. The 
previous two strong cases of reverse shock emission (GRB 
990123, Akerlof et al. 1999; and GRB 021211, Fox et al. 
2003, Li et al. 2003) all belong to Type II. Visual inspec- 
tion of the early IR lightcurve of GRB 021219a (Blake et 
al. 2005, see also Fig|2J) indicates that it is a clear Type 
I case. Below we will model the lightcurve in detail and 
show that the data are indeed consistent with such an 
explanation. 

Following the standard afterglow model for a fire- 
ball interacting with a constant density medium (e.g., 
Piran 1999), we write down the cooling frequency v , 
the typical synchrotron frequency and the maxi- 
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83mJy £^ so ,54%' j _2"'0 
7(p — 2)/[2(p — 1)], E iso is the isotropic energy of the 
outflow, e e and £b are the fractions of the shock energy 
given to electrons and to magnetic fields in the forward 
shock, respectively, n is the number density of the ex- 
ternal medium, p ~ 2.4 is the power-law distribution 
index of the shocked electrons. Hereafter t = i bs/ (1 + z) 
denotes the observer's time corrected for the cosmolog- 
ical time dilation effect, and t& is in unit of day. The 
superscript "f" ("r") represent the FS (RS) emission, re- 
spectively. We assume that e c and the electron spectral 
index p are essentially the same for both the FS and RS 4 , 
but we will allow different £b values for both regions. 
One reason for this assumption is that the magnetic field 
generated in the internal shock phase may have not been 
dissipated effectively in a short time, and would play a 
dominant role in the reverse shock region (Fan et al. 2004 
and the references listed therein). In this Letter, we will 
denote £ B = £b and e B = 1Z B £ B, where TZb is the ratio 
of the magnetic field in the RS emission region to that 
in the FS emission region (Zhang et al. 2003). Previ- 
ous analyses indicate that at least for some bursts (e.g. 
GRB990123 and GRB021211) the RS emission region is 
more magnetized than the FS region (e.g. Fan et al. 
2002; Zhang et al. 2003; Kumar & Panaitescu 2003). 

As shown in Fig[5](data taken from Blake et al. (2005), 
only the richest i^ s -band data are plotted), the time 
when the RS crosses the ejecta ((l+z)t x ) is about 30 
minutes after the trigger, which is much longer than 
T90 ~ 520s. So, the RS is non-relativistic (e.g., Sari 
& Piran 1995; Kobayashi 2000; Kobayashi & Zhang 
2003). Further evidence for a non-relativistic RS is the 
rapid increase of the early afterglow lightcurve (e.g., 
Kobayashi 2000; Kobayashi & Zhang 2003). 734.x « 

4 If they are different, additional corrections are needed, see 
Zhang et al. (2003) and Zhang & Kobayashi (2005) for details. 
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Fig. 2.— Modeling the K s -band afterglow data of GRB 041219a. 
The data are taken from Blake et al. (2005). The earliest data are 
prompt emission, which are not included in our fit. The dash- 
dotted line is the theoretical lightcurve. Both the reverse shock 
and the forward shock emission components are included. The 
best fit parameters arc marked in the text box. 



(?y/r x + T x /rj)/2, the LF of the decelerated outflow rel- 
ative to the initial LF at t x , can be estimated by solving 

cdt/dA w (l-/3r 3 ){l-r ? /[r3(4734+3)]}/(/3 7) -/3r 3 ) (e.g., 
Sari & Piran 1995; Fan et al. 2004) numerically. Here T 3 
is the bulk LF of the shocked ejecta, rj is the initial bulk 
LF of the outflow, T x is the LF of the decelerated ejecta 
at t x , Pa is the corresponding velocity (in unit of c) of 
the LF T A . For GRB 041219a, one has t x ~ 4T 90 /(l + z), 
so that 734.x — 1 is much smaller than 1, i.e., the RS is 
non-relativistic. In such a case, t x can be approximated 



as 



1/3 -1/3 -8/3 
iso,54 n 



64s £,„'„-,//,, ' -772.5 



(1) 



The typical frequency of the RS emission is 

i4(t x ) =ft B (734,x - i) 2 *4(tx)/(r x - 1) 2 , (2) 

Following Zhang et al. (2003), we have 

u\ « ll^vl K,mJt*) « rfRsFt, m»(*x). (3) 
Generally, the K s band flux satisfies F VK (t x ) « 

^ m ax(^)[^K s /^(tx)]- (p - 1)/2 . 

Assuming z = 1, p = 2.4 and -Eiso. 54 = 1 (consistent 
with the gamma-ray fluence and our assumed redshift), 
e , n and 77 can be constrained by the following two con- 
ditions. (1) At td ~ 0.14, v { m ~ j/k 3 and F^ max ~ 0.6 
mJy; (2) t x is about 30/(l+z) minutes. We then get 

£ C -0.2£ B V 2, n - 5 x 10" 5 £ B 1 _ 2 , 7y~ 380e^-2- ( 4 ) 

It is interesting to see that since £ e < 1, we have £b > 
1.6 x 10~ 5 , n < 3 x 10~ 2 cm~ 3 and 77 > 170. The param- 
eter TZb can be constrained by noticing F„ Ka [(1 + z)t x ] ~ 
2.4 mJy and by taking 77 ~ 380, which reads 



n B ~ 3[( 7: 



34, X 



l)/0.2] 



-14/17 



(5) 



which hints that the reverse shock region is mildly mag- 
netized. 

Following Fan et al. (2005) , the forward-reverse shock 
emission has been calculated numerically. Since the RS 
is non-relativistic, the spreading of the ejecta (e.g. Piran 
1999) has been taken into account. The fits to the K s 
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band data are presented in Fig. [21 It is found that the 
data can be well modeled with the following parameters: 
z = 1, E isoM = 1, e e = 0.45, e B = 0.01, K B = 2.8, 
r\ = 300, n = 10~ 4 cm -3 . These are consistent with the 
analytical estimates above. 

There are three radio data points available (Soderberg 
& Frail 2004; van der Horst et al. 2004a, 2004b). The 
8.5GHz flux at 1.1 day is 0.45mJy, and the 4.9 GHz 
flux at 1.6 day and 2.6 day are 0.2mJy and 0.34mJy, 
respectively. By taking our best fitted parameters, the 
corresponding FS fluxes are ~ (0.04, 0.04, 0.05)mJy, 
respectively, too low to interpret the data. When we 
consider the RS contribution to the radio band (e.g. 
Sari & Piran 1999; Kobayashi & Zhang 2003; Gou et 
al. 2004), the over all corresponding fluxes become 
~ (0.60, 0.37, 0.18)mJy, respectively, roughly match- 
ing the data. Therefore, the radio data also support the 
reverse-forward shock interpretation. 

4. SUMMARY & DISCUSSION 

The prompt optical/IR observations of GRB 041219a 
(Vestrand et al. 2005; Blake et al. 2005) offer a great op- 
portunity to diagnose the unknown GRB ejecta compo- 
sition. We have shown that the prompt optical emission 
tracking the gamma-ray emission profile may be consis- 
tent with the picture that the ejecta is neutron rich, and 
that the optical emission is powered by the proton shells 
interacting with the neutron decay products at a distance 
farther away from the central engine than the typical in- 
ternal shock radius. 

By modeling the K s band early afterglow lightcurve, 
we identify a reverse shock emission component, which is 
clearly separated from the forward shock emission com- 



ponent. Such a rebrightening (Type I) lightcurve has 
been expected by Zhang et al. (2003) to be more com- 
mon if the RS is not strongly magnetized. Indeed, de- 
tailed modeling indicates that TZb in GRB 041219a is 
at most mild, in contrast with GRB 990123 and GRB 
021211 (e.g., Fan et al. 2002; Zhang et al. 2003; Ku- 
mar & Panaitescu 2003). This is also consistent with 
the neutron-rich picture conjectured in interpreting the 
prompt optical emission. The mild magnetization of the 
ejecta may be due to magnetic field generation during 
the internal shock phase. The radio data also support 
our reverse-forward shock interpretation. 

It is interesting to note that for the three bursts with 
reverse shock identification (GRB 990123, GRB 021211 
and GRB 041219a), the inferred number densities of in- 
terstellar medium are typically lower than the standard 
value for ISM n ~ 1 cm" 3 (see also Kumar 2004) — For 
GRB 990123, n ~ 10~ 3 cm" 3 (e.g., Wang et al. 2000; 
Panaitescu & Kumar 2002; Nakar & Piran 2004); For 
GRB 021211, n~ lO^-lO" 2 cm~ 3 (e.g., Wei 2003; Ku- 
mar & Panaitescu 2003); For the current GRB 041219a, 
n ~ 10~ 4 cm~ 3 . The reason is unclear, but certain se- 
lection effects may play a role. 
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